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Abstract 

Here we use polynomial chaos framework to design controllers for linear parameter varying (LPV) dynamical systems. We 
assume the scheduling variable to be random and use polynomial chaos approach to synthesize the controller for the resulting linear 
stochastic dynamical system. The stability of the LPV system is formulated as an exponential mean-square (EMS) stability problem. 
Two algorithms are presented that guarantee EMS stability of the stochastic system and correspond to parameter dependent and 
independent Lyapunov functions, respectively. LPV controllers from the polynomial chaos based framework is shown to outperform 
LPV controller from classical design for an example nonlinear system. 


I. Introduction 

Linear parameter varying (LPV) systems are of the form 

X = A{p)x + B{p)u, (1) 

where system matrices depend on unknown parameter p{t), which is measurable in real-time HI, ||2l. Many nonlinear systems 
can be transformed to LPV systems and control systems can be designed using parameter dependent convex optimization 
problems. Typically, parameter dependent quantities are approximated using a known class of functions such as multilinear 
basis functions of p, linear fractional transformations of system matrices, or by gridding the parameter space. Both these 
approaches result in solution of a finite, but possible large, number of linear matrix inequalities (LMIs). Further, the choice of 
the basis functions or the resolution of the grid could lead to conservatisms in the design. Clearly, there is a tradeoff between 
problem size and conservatism in the design 0. 

Fujisaki et al. 11 addressed the computational complexity of such problems by presenting a probabilistic approach to solve 
these problems, via a sequential randomized algorithm, which significantly reduces the computational complexity. Here the 
parameter p{t) is assumed to be bounded i.e. p{t) S Pp C and is treated as a random variable, with a distribution fp{p) 
defined over Vp. The LPV synthesis problem is solved by sampling Vp and solving the sampled LMIs using a sequential- 
gradient method. As with any probabilistic algorithm, there is a tradeoff between sample complexity and confidence in the 
solution. Often, a large number of samples are required to generate a solution with high confidence. Also, the LMIs depend 
only on p{t) and not in p{t) as it is in classical LPV formulation. 

This paper is motivated by the work of Fujisaki et al. and is based on the idea of treating p as a random variable. Therefore, 
by substituting p = A in the system equation, we get 

X = A{A.)x + B{A.)u, (2) 

where A G is a vector of uncertain parameters, with joint probability density function /a(A). Matrices A.(A) G 
B (A) G are system matrices that depend on A. Consequently, the solution x := x{t,A.) G K" also depends on A. 

Like in ||4l we ignore temporal variation in the parameter and thus treat A as random variables. Thus, we now study the 
system in ([T]) as a linear time invariant system with probabilistic system parameters. The LPV control design objective is then 
equivalent to designing a state-feedback law of the form u = K{A.)x, which stabilizes the system in some suitable sense, 
where K{A) G Thus, we are looking to obtain a parameter dependent gain K{A) that stabilizes the system in (|2]). 

The closed-loop system is then 

X = [A{A) + B{A)K{A)]x, 

= Acix. (3) 

There are two distinct differences between the work presented here and that in ||4|. We do not use a randomized approach to 
solve the stochastic problem, and thus don’t have issues related to confidence in the solution. In our approach, the stochastic 
problem is solved using polynomial chaos theory, which is a deterministic approach as described later. In addition, stability 
of the LPV system is formulated as an exponential mean square stability problem for the corresponding stochastic system. In 
a, stability of the LPV system is formulated in the probabilistic sense. Computationally, the polynomial chaos framework is 
superior to sampling based approach in propagating uncertainty a, and hence we can expect a computational advantage in 
using this framework to solve the stochastic formulation. 
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Main contributions of this paper are two LPV controller synthesis algorithms with parameter dependent and independent 
Lyapunov functions respectively. They are presented as theorem 1 and 2. The paper is organized as follows. We hrst provide 
a brief background on polynomial chaos theory and show how it is applied to study linear dynamical systems with random 
parameters. This is followed by conditions for exponential mean-square stability in the polynomial chaos framework for 
closed-loop systems with parameter dependent controller. This leads to theorem 1 and 2. The paper ends with an example that 
highlights the superiority of the polynomial chaos approach over the classical LPV design approach. 

11. Polynomial Chaos Theory 

Polynomial chaos (PC) is a non-sampling based method to determine evolution of uncertainty in dynamical system, when 
there is probabilistic uncertainty in the system parameters. Polynomial chaos was hrst introduced by Wiener 0 where Hermite 
polynomials were used to model stochastic processes with Gaussian random variables. It can be thought of as an extension 
of Volterra’s theory of nonlinear functionals for stochastic systems Q, 0. According to Cameron and Martin 0 such an 
expansion converges in the Li sense for any arbitrary stochastic process with hnite second moment. This applies to most 
physical systems. Xiu et al. Eol generalized the result of Cameron-Martin to various continuous and discrete distributions 
using orthogonal polynomials from the so called Askey-scheme ini and demonstrated Ci convergence in the corresponding 
Hilbert functional space. The PC framework has been applied to applications including stochastic huid dynamics Ha-Ill, 
stochastic hnite elements 0, and solid mechanics na, m, feedback control Ha-GQl and estimation IHI. It has been shown 
that PC based methods are computationally far superior than Monte-Carlo based methods 0, Ha, na-Hl- see Ea for 
several benchmark problems. 

A general second order process X{u!) € P) can be expressed by polynomial chaos as 

OO 

X(a;) = ^x*0i(A(a;)), (4) 

where uj is the random event and (j)i{A.{u!)) denotes the polynomial chaos basis of degree p in terms of the random variables 
A(a;). is a probability space, where il is the sample space, P is the cr-algebra of the subsets of H, and P is the 

probability measure. According to Cameron and Martin 0 such an expansion converges in the Li sense for any arbitrary 
stochastic process with hnite second moment. In practice, the inhnite series is truncated and X{u!) is approximated by 

N 

X{uj) ^ X{uj) = y^^Xi(l)i{A{uj)). 

i=0 

The functions are a family of orthogonal basis in Li{ft,P, P) satisfying the relation 

E := [ <^,(A)</),(A)/a(A) dA = h^Ay-, (5) 

where A^ is the Kronecker delta, hi is a constant term corresponding to (j)^f^{A) dA, is the domain of the random 
variable A(a;), and /a(A) is a probability density function for A. Table shows the family of basis functions for random 
variables with common distributions. 


Random Variable A 

0i(A) of the Wiener-Askey Scheme 

Gaussian 

Hermite 

Uniform 

Legendre 

Gamma 

Laguerre 

Beta 

Jacobi 


TABLE I: Correspondence between choice of polynomials and given distribution of A(u;) □a. 


A. Application to Dynamical Systems with Random Parameters 

With respect to the dynamical system dehned in Q, the solution can be approximated by the polynomial chaos expansion as 

N 

x{t,A) Ki x{t,A) = '^Xi{t)cj),{A), (6) 

i^O 

where the polynomial chaos coefficients Xi € Define ^(A) to be 

^ = ^(A) := ((/>o(A) (/)Ar(A))^ , and 

$„ = #„(A) :=#(A)(g)I„, 


(7) 

( 8 ) 





where J„ S is identity matrix. Also define matrix X G with polynomial chaos coefficients Xi, as 

X = [ato • • • xn ] ■ 


This lets us define x{t, A.) as 


x{t,A) :=X(f)$(A). 

Noting that x = vec (x), we obtain an alternate form for (|9ll, 

X = vec (x) = vec (A$) = vec (I„X#) = ® J„)vec (X) = ^^Xpc, 

where Xpc := vec (X), and vec (•) is the vectorization operator ll23l . 

Since x from (fTOl i is an approximation, substituting it in (O we get equation error e, which is given by 

G . — X A(‘i(^A')x — ^^Xpc A.c/(tCpc- 
Best £2 approximation is obtained by setting 

(ec^i) := E [e4ii] =0, for * = 0,1, • • • , 


(9) 

( 10 ) 


( 11 ) 

( 12 ) 


* 


Xpc — E 


Xpc — E 


^nAcl^ 

-1 


* 




or Xpc — ApcXpc- 


where and Ad depend on A as defined earlier. 

We will need the following result in the rest of the paper. 

Proposition 1: For any vector v G and matrix M G 

M(v'^ 0 In) = (v'^ (g) Im)(lN+l ^ M), 


where J* is identity matrix with indicated dimension. 
Proof: 


M(u^ ® In) = (1 0 M){v^ ® In) 

= v^®M= ( u '^/ tv + i ) ® (ImM) 

= (g) 7„)(7Ar+i 0 M). 


(13) 

(14) 


(15) 


III. Controller Synthesis 

The controller gain K{A) can be introduced in the polynomial chaos framework by substituting Ad := A{A) + B{A)K{A), 
in (O to get 


^pc - ^ 


* 


^nA^^n 


^nBK^i 


Xri 


(16) 


Polynomial chaos expansion of K{A) can be written as 


N 


k{A) = Y,km^),k, g 


2 = 0 


— [00-^r] 


(j)Nlr) 


Kn 


K 


N 


= {^^<E>Im)VK = ^lVK, 

where Vk € is the vertical stacking of Ki. The expression BK^^ in (fTST l can be simplified using (fTST i 


(17) 




where Vk '■= In+i ® ^k- Therefore, 


^pc - E 


* 






VjC ) Xpc. 


(18) 
























Recall that for the dynamical system in Q, the equilibrium solution is said to possess exponential stability of the m* mean 

if 3 A > 0 and constants a > 0,/3 > 0 such that ||a:o|| < A implies Vf > to liSjl . ll25l 

E[||a;(f;a:o,fo)C]</?E[||xo||™]e-“(‘-“>). (19) 

It can be shown ll26l that the dynamical system in Q, with random variables A, is exponentially stable in the 2"‘* mean, or 

exponentially stable in the mean square sense (EMS-stable), if 3 a Lyapunov function V{x) := x^Px, with P = > 0, 

and a > 0 such that 


V 


< —aE [V]. 


( 20 ) 


Theorem 1: The closed-loop system (|3ll is EMS-stable with controller K{A.) if 3 P = P^> 0 and a > 0 such that 


s> 4*^ 

T 

+ E 

^uA^l 

y + w^E 


+ 

m 


W-f a^E 



< 0 , 


( 21 ) 


where W := In+i (g)W,y:= In+i (SiY,Y := p-\ and W := VkY. 

Proof: With V{x) := x^Px, and P = P^ > 0, E = x^Px -f x^Px. The term x^Px can be approximated by the 
polynomial chaos expansion as 


x^Px 


pc 






= xfJ<&^A^ P<^i +Vl<^ 


pc 




P^i 


n I ^pc, 


Using (fTSl l we can write P$^ = where V := Tat+i <8) P. Substituting them in x^Px we get 


x^ Px « xf^ ( 


pc 


K^m{N+l)^m 




The Lyapunov function V := x^Px can be written as 

V := x^ 

Therefore, E V < — aE [U] is equivalent to 


V := X^PX = Xp^^riP^n^pc = Xp^^n^'^PXpc. 


( 22 ) 


$ 




P + V^E 


^m(Ar-t-l) 




P + (*)^ ^ —aE 




P, 


where (*)^ are the symmetric terms. The above BMI can be convexified using the well known substitutions lIZTl Y := P ^ 
and W := VkY. These substitutions can be written in terms of V,VK,y, and W as 


W = Ttv+i ® = In+iIn+i ® VkY 

= (In+i ® Vk){IN+ i ® Y) 

= VKy. 

It is also straightforward to show P = y~^ and Vk = Substituting these in the above BMI, and pre-post multiplying 

by y, we get the result. ■ 

Theorem 2: The closed-loop system Q is EMS-stable with controller K{A) if 3 P(A) = P^(A) > 0 and a > 0 such that 

3- MiE + + MaW -f ay Mo < 0, (23) 

where 


P(A) :=#^(A)P#„(A), 

P •= In+1 ^ Pqi Pq = Pq > 0) 


Mo = E 
Ml = E 
Ma = E 


* AS>^ 

^n^n^rL-^^n i 




























Proof: Define Lyapunov function V{x) := x'^P{A)x, with -P(A) := $^(A)P$„(A), V := In+i ® Po, and Pq = 
Pq > 0. The form for -P(A) is motivated by the literature on sum-of-square representation of matrix polynomials 1281 . l2^ . 
which ensures P(A) > 0. The Lyapunov function can be simplified as 

V(x) = x’^P{A.)x 

= Xp^^n^l (In+I < 8 ) Po)^n^lxpc 
= ® In){lN+l ® ® In)Xpc 

= ® ® Po)Xpc 

0 P o)'^pc 

V = jFPx + x^Px. The term xFPx can be approximated by the polynomial chaos expansion as 


x'^P(A)x « xl^ 

We next show that 

V^n^n = {In+1 0 -Po)(^ 0 0 In) 

= (W®Po)((^^^)0-fn) 

= 

= Jat+i (g) InPo 

= ((^$^) (g) Ir^In){IN+l 0 Po) 

= 


Therefore, 


xFP{A)x 

= Xpc tCpc- 


Therefore, E 


V 


< —aE [V] is equivalent to 


MfP + PMi + VlcM'^r + VM 2 VK + aMoV < 0, 


which can be convexified as in Theorem 1 to obtain the result. 


IV. Example 


Here we consider the control of the following nonlinear system 



0 

-1 


1 

(1-xf) 




(24) 


The above systems is the Van der Pol oscillator with a control input. We transform it to an LPV system by introducing the 
parameter p := 1 — x^. The objective is to design a state feedback controller K{p) that will quadratically stabilize the above 
system. We restrict stabilization of the set defined by at G [—5,5]^. Therefore, p € [—24,1]. For the PC LPV algorithm, we 
assume p = A G Z^[_ 24 .i]. ^ uniformly distributed random variable over [—24,1]. Fig.([T]i shows the state and control trajectories 
of (l24l i with three control systems /Clxl -f^LPv(p) and /CpcLPv(A), and were designed with a = 1 in the following manner: 
• iCuri, from linearized dynamics (Alti,Plti) about (0,0), satisfying 


Y ltiA^j + AltiT^lti + + .BltiWlti + oiY lti < 0. 


• iCLPv(p), from LPV dynamics 



1 

P 


,Blpv(p) := 


0 

1 


Alpv(p) : 












satisfying 

y LPv(Pfe)^LPv(Pfc) + ^LPv(Pfe)^ LPv(Pfe) + ^LPv(Pfe)-®LPv(Pfe) + -BlPV (pfe) W'lpV (Pfe) + OiY LPv(pfe) < 0, 
where 

l^LPv(pfe) := 1^0 + Pfc>"i > 0, y, = Y[, 

and pk are the samples from W[_ 24 ,i]- 
« i^pcLPv, from theorem 2, assuming p = A G ^_ 24 ,i]- 

Fig.([T]i shows the state and control trajectories of the nonlinear closed-loop system for initial condition (5, 5). i^pcLPv is 
designed with first order polynomial chaos expansion and several -f^LpyS are designed with 2, 5,10 and 50 samples from 

^- 24 , 1 ]- 

For this problem, we make the following observations. 

1) Increasing the order of the polynomial chaos expansion does not significantly improve controller performance. We are 
able to achieve high performance with very low order polynomial chaos expansion. 

2) As seen from fig.© increasing the number of samples in the design of FFlpvj improves the performance, but doesn’t 
quite reach the performance of i^pcLPV- The computational times for iTupv synthesis are as follows; 


Controller 

Synthesis Time (s) 

JVpcLPV (first order PC) 

0.3556 

IClpv (2 samples) 

0.4008 

-f^LPV (5 samples) 

0.5463 

TFlpv (10 samples) 

0.6723 

TFlpv (50 samples) 

1.9943 


Thus FFpcLPV has a clear advantage over sampled based FFlpv design in terms of controller performance and computational 
complexity. 

3) Both iTpcLPv and FFlpv outperform JFlti as expected. 

The controllers were synthesized in MAT LAB using CVX ED- 

V. Summary 

In this paper we presented a new framework to synthesize LPV controllers using polynomial chaos framework. This framework 
builds on the probabilistic representation of the scheduling variable and the synthesis was done by treating the LPV system as 
a stochastic linear system. Two synthesis algorithms were presented which correspond to parameter dependent and independent 
Lyapunov functions. The algorithms were tested on a nonlinear dynamical system and outperformed controllers synthesized 
using classical LPV design techniques. 
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